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ABSTRACT 

Black hole accretion and jet production are areas of intensive study in astrophysics. Recent work has found 
a relation between radio luminosity, X-ray luminosity, and black hole mass. With the assumption that radio 
and X-ray luminosity are suitable proxies for jet power and accretion power, respectively, a broad fundamen- 
tal connection between accretion and jet production is implied. In an effort to refine these links and enhance 
their power, we have explored the above relations exclusively among black holes with direct, dynamical mass- 
measurements. This approach not only eliminates systematic errors incurred through the use of secondary mass 
measurements, but also effectively restricts the range of distances considered to a volume-limited sample. Fur- 
ther, we have exclusively used archival data from the Chandra X-ray Observatory to best isolate nuclear sources. 
We find logL R = (4.80±0.24)+ (0.78 ±0.27)logM BH + (0.67 ± 0.12)logL x , in broad agreement with prior 
efforts. Owing to the nature of our sample, the plane can be turned into an effective mass predictor. When 
the full sample is considered, masses are predicted less accurately than with the well-known M-a relation. If 
obscured AGN are excluded, the plane is potentially a better predictor than other scaling measures. 
Subject headings: black hole physics — galaxies: general — galaxies: nuclei — galaxies: statistics 



1. INTRODUCTION 

Accretion onto black holes has many observable conse- 
quences, including the production of relativistic jets. The phe- 
nomenon of jet production appears to be universal, as such 
jets are observed both in Active Galactic Nuclei (AGN) and 
stellar-mass black hole systems as well as in neutron stars, 
white dwarfs, and even young stellar objects. For black hole 
sources, the length scales and relevant timescales of jets ap- 
pear to approximately scale with mass over 8 orders of mag- 
nitude, giving rise to the possibility that jet production mech- 
anisms scale with mass, similar to the way that accretion disk 
properties scale. The mechanism by which jets are driven 
from black holes, however, remains observationally elusive. It 
remains one of the most compelling and important problems 
in astrophysics, particularly in high energy astrophysics. The 
impact of relativistic jets on the interstellar medium (Gallo 
et al. 2005b), and large-scale structure in clusters of galaxies 
(Allen et al. 2006; Fabian et al. 2003; McNamara et al. 2006), 
has become dramatically clear in the era of imaging and spec- 
troscopy with Chandra. 

Virtually all theories of jet production tie the jet to the ac- 
cretion disk directly or indirectly (see, e.g., Lynden-Bell 1978; 
Blandford & Payne 1982, see also van^Putten 2009). Thus, 
there is a broad expectation that jet properties might depend 
on the mass accretion rate (M) through the disk. The black 
hole spin parameter (a = cJ/GM 2 ; < a < 1) may also be 
an important factor if the black hole and accretion disk are 
linked through magnetic fields (Blandford & Znajek 1977). 
The spin is also important for accretion disk jet-launching be- 
cause the inner radius of the accretion will decrease, thus in- 
creasing the launch velocity. This idea may find some support 
in the dichotomy between radio-loud and radio-quiet AGN 
(Sikora et al. 2007). The high flux of stellar-mass black holes 
facilitates spin constraints with current X-ray observatories; in 
those systems, the most relativistic jets appear to be launched 
by black holes with high spin parameters (Miller et al. 2009). 

4 Chandra Fellow 



One means by which jet production can be examined is 
to explore correlations between proxies for mass inflow and 
jet outflow. In stellar-mass black holes, it was found that 
radio emission and X-ray emission are related by Lr °c L x 7 
(Gallo et al. 2003). This correlation was quickly extended to 
also include super-massive black holes in AGN, resulting in 
the discovery of a "fundamental plane" of black hole activ- 
ity (Merloni et al. 2003; Falcke et al. 2004, also see Merloni 
et al. 2006). The plane can be described by logLfl = 7.33 + 
0.601ogLx +0.781ogM B H, with a scatter of o R = 0.88 dex 
(where Lr is v = 5 GHz nuclear radio luminosity in units 
of ergs -1 , Lx is E = 2-10 keV nuclear X-ray luminosity in 
units of ergs -1 , and Mbh is the black hole's mass in units 
of M . : Merloni et al. 2003; Falcke et al. 2004). Several re- 
cent works have revisited the original findings with slightly 
different focuses. Kording et al. (2006) found that sources 
emitting far under their Eddington limits followed the relation 
more tightly. Wang et al. (2006) found differences in the rela- 
tionship for radio-loud and radio-quiet AGNs. Li et al. (2008) 
used a large sample of SDSS-identified broad-line AGNs to 
study a similar relation at lower-frequency (1.4 GHz) radio 
luminosity and softer-band (0.1-2.4 keV) X-ray luminosities. 
Yuan et al. (2009) limited the sample to those sources with 
Lx/^Edd < 10 6 based on predictions that the correlation be- 
tween radio and X-ray luminosity steepens to Lr °< L x 2i at 
low accretion rates (Yuan & Cui 2005). 

It is difficult to overstate the potential importance of the 
fundamental plane; it suggests that black holes regulate their 
radiative and mechanical luminosity in the same way at any 
given accretion rate scaled to Eddington, m = M/Medd- In the 
context of models that assume jet properties do scale simply 
with mass (e.g., Falcke & Biermann 1995; Heinz & Sunyaev 
2003), the fundamental plane can even be used to constrain 
the nature of the accretion inflow. At present, radiatively- 
inefficient inflow models for X-ray emission, and models as- 
sociating X-ray flux with synchrotron emission near the base 
of a jet, are both consistent with the fundamental plane. 

To use the fundamental plane as a tool and a diagnostic 
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instead of as an empirical correlation, however, it must be 
sharpened. Black hole masses represent a significant source 
of uncertainty and scatter in the fundamental plane (Merloni 
et al. 2003; Kording et al. 2006). In this work, we have 
constructed a fundamental plane using only black holes with 
masses that have been dynamically determined, the so-called 
M-o black holes (see Giiltekin et al. 2009b). Unlike prior 
treatments, our X-ray data is taken from a single observa- 
tory and predominately from a single observing mode, and we 
have conducted our own consistent analysis of the data. We 
analyzed every archival Chandra X-ray observation of black 
holes with a dynamically-determined mass. Radio data were 
taken from archival observations reported in the literature. By 
using a sample of black holes with dynamical masses, we may 
probe the fundamental plane without subjecting the analysis to 
the systematic errors inherent in substituting scaling-relation- 
derived quantities for black hole masses. 

In Section 2 we describe the sample of black holes used in 
this work. We detail our X-ray data reduction and spectral fits 
in Section 3. Our fitting methods and results are presented in 
Section 4. We discuss our results in Section 5 and summarize 
in Section 6. 

2. SAMPLE OF BLACK HOLES 

We get Mbh from the list of black hole masses compiled 
in Giiltekin et al. (2009b), adopting the same distances as 
well. This sample of black hole masses includes measure- 
ments based on high spatial resolution line-of-sight stellar 
velocity measurements (e.g., Giiltekin et al. 2009a), stellar 
proper motions in our Galaxy (Ghez et al. 2008; Gillessen 
et al. 2008), gas dynamical measurements (e.g., Barth et al. 
2001), and maser measurements (e.g., Miyoshi et al. 1995). It 
does not include reverberation mapping measurements, which 
are direct measurements of mass but are secondary in that they 
are normalized to the other measurements via the M-a rela- 
tion (e.g., Peterson et al. 2004; Onken et al. 2004). From those 
available black hole masses we use only the measured black 
hole masses used in their M-a fits — not upper limits and not 
the "omitted sample," which contains a list of masses with po- 
tential problems because (1) masses were listed as tentative 
by the original study, (2) there was no quantitative analysis 
of how well the original study's model fit the data, or (3) the 
quantitative analysis of the goodness of fit was poor. We re- 
duce available Chandra data and present X-ray luminosities 
for this collection of potentially problematic masses, but we 
do not use them in our fits. Thus, we use only the black hole 
masses with the most reliable measurements. 

One of the benefits of using this sample is that most of the 
distances to the galaxies are less than 30 Mpc. This distance is 
close enough that interestingly low X-ray and radio luminosi- 
ties will still be measurable. So while this is not a true volume- 
limited sample, it is insensitive to the potential biases arising 
from, e.g., a sample limited by X-ray flux. Unlike other sam- 
ples, however, our sample may be biased to very low nuclear 
luminosities. The contamination from a bright AGN typically 
causes problems in determining stellar mass-to-light ratio at 
the center so that most galaxies selected for dynamical mea- 
surement do not contain bright AGN. Spiral galaxies, which 
are less massive on average than early-type galaxies, may also 
be underrepresented in this sample, and thus low-mass black 
holes may also be underrepresented. 

The radio data we use are 5 GHz peak power measurements 
from the Ho (2002) compilation of nuclear radio sources. The 
data were compiled to probe whether there was a correlation 
between Mbh and L^and thus are ideal for our purposes. 

3. X-RAY ANALYSIS 



3.1. X-ray Data Reduction 

The high spatial resolution of Chandra enables nuclear 
emission to be isolated best compared to other X-ray obser- 
vatories. We used Chandra archival data to obtain accurate 
measurements or tight upper limits of the flux between 2 and 
10 keV for all galaxies in our sample. This energy range was 
chosen to probe accretion power rather than total power in- 
cluding any contaminating diffuse emission and for ease in 
comparison to previous fundamental plane work. 

For each source we used a circular extraction region posi- 
tioned at the brightest point source that was consistent with 
the center of galaxy determined by 2MASS images. There are 
no nuclear point sources in NGC 1399 and NGC 4261, which 
we handled slightly differently as described below. For ex- 
traction of background spectra, we typically used an annular 
region with inner radius slightly larger than the source region 
radius. The outer radius was made large enough to encompass 
a significant number of counts. Point sources were excluded 
from the background region. When there were a large number 
of point sources in the annular region surrounding the source, 
a different region was used, usually an off-nuclear circle. In 
these cases we selected a region where the background looked 
to be similar to that surrounding the source. 

Two cases require special attention to contamination 
from non-nuclear X-ray emission: NGC 0224 (M31) and 
NGC 4486 (M87). NGC 0224 has two bright point sources 
near the center of the galaxy, but neither is the galaxy's central 
black hole, from which the emission is too dim to be detected 
above the background to high significance (Lx < 10 36 erg s _1 
at at assumed distance of D = 0.8 Mpc Garcia et al. 2005; Li 
et al. 2009). NGC 4486 is well known for its prominent jet 
with several knots. These knots are apparent in the Chandra 
images. In the high spatial resolution Hubble Space Telescope 
(HST) images, one knot is very close (0"85 Harris et al. 2006) 
to the central engine. As the knot has grown brighter in the 
optical by a factor of ~ 100 over the last ~ 10 years, mea- 
surements of the core X-ray flux become increasingly con- 
taminated by the knot. We chose the archival Chandra data 
set where the knot was most readily distinguishable from the 
core. 

For the galaxies NGC 1399 and NGC 4261, there is no dis- 
cernible point source at their nuclei, which are dominated in 
X-rays by hot gas. For these two sources, we attempt to mea- 
sure a hypothetical point source at the center. We use a circu- 
lar region at the center of the diffuse X-ray emission for source 
extraction with an annular background extraction region im- 
mediately adjacent. For both these sources, X-ray point source 
flux could not be inferred above the background, and they are 
listed as upper limits in Table 3. 

Data reduction followed the standard pipeline, using the 
most recent Chandra data reduction software package (CIAO 
version 4.1.1) and calibration databases (CALDB version 
4.1.2). Point-source spectra were extracted using the CIAO 
tool psextract. Because all observations of interest were done 
with the Advanced CCD Imaging Spectrometer (ACIS), we 
ran psextract with the mkacisrmf tool to create the response 
matrix file (RMF) and with mkarf set for ACIS ancillary re- 
sponse file (ARF) creation. 

3.2. X-ray Spectral Fitting 

We modeled the reduced spectra using XSPEC12 (Arnaud 
1996). If binning the spectra in energy so that each bin con- 
tained a minimum of 20 counts resulted in five or more bins, 
we did so and used % 2 statistics; otherwise we did not bin 
the data and used C-stat statistics (Cash 1979). Each spec- 
trum was modeled with a photoabsorbed power-law model. 
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If such a model did not adequately fit the spectrum for data 
sets that were strong enough to support a more complicated 
model, we added additional model components. Galaxies 
that were identified as Seyfert 2 or transitional Seyferts in 
Veron-Cetty & Veron (2006) were modeled with a partially 
photoabsorbed power-law, representing intrinsic absorption 
plus another photoabsorbed component, representing Galac- 
tic absorption. Galaxies with obvious diffuse hot gas towards 
their nucleus were modeled with photoabsorbed Astrophys- 
ical Plasma Emission Code (APEC Smith et al. 2001) and 
power-law components. Regardless of the continuum model, 
for spectra that showed an obvious Fe Koc line, we added a 
Gaussian for each line. All spectra were fit from E = 0.5 to 
10 keV. 

We considered a model successful if it yielded a reduced % 2 
of % 2 / v < 2 and if the spectrum between E = 2 and 10 keV 
was adequately described. The total flux between E = 2 
and 10 keV, Fx,t t, was determined from the model and the 
la errors derived from covariance of the model parameters. 
We then calculated the unabsorbed flux arising from just the 
power-law component between E = 2 and 10 keV, Fx. That 
is, we de-absorbed the flux and removed contributions from 
lines and other model components. We assume the fractional 
error in Fx, tot is the same as in Fx- 

For sources that did not constrain the flux from the cen- 
tral point source, we used the total count rate between E = 
0.5 and lOkeV to calculate the 3a (99.7% confidence) up- 
per limit to Fx with PIMMS assuming a power-law with in- 
dex r = 2 and with Galactic absorption determined from the 
Leiden/ Argentine/Bonn survey of Galactic HI (Kalberla et al. 
2005; Hartmann & Burton 1997; Bajaja et al. 2005) using the 

HEASOFT ftool "N H "- 

Because we are ultimately interested in an accurate mea- 
surement of Fx, it is more important that our models charac- 
terize the spectrum well over the 2 to 10 keV band than it is 
to reproduce the underlying physics. We tested this approach 
by fitting several different models to the same spectrum and 
recovered consistent values for Fx. The results of fits are dis- 
played in Table 2, and we show four example spectra with 
models in Figure 1 . 

For many galaxies, multiple Chandra observations were 
available in the archive. We reduced and analyzed the avail- 
able data and censored the resulting data by (1) choosing those 
that yielded flux detections as opposed to upper limits, (2) 
choosing those with smaller values of % 2 /v, (3) preferring 
higher precision measurements over lower precision, and (4) 
observed more closely in time with the available radio data 
since variable sources will have Lr and Lx change in concert 
on the fundamental plane (see Merloni et al. 2006). 

We compare our results with results from the literature for 
the same data sets in Figure 2. The literature values were 
scaled to our assumed distances and, in some cases, converted 
to the 2-10 keV band with PIMMS and the published spectral 
fits. The comparison reveals good agreement with no partic- 
ular bias with exception of a single outlier, NGC 1068. We 
expand on NGC 1068 and Compton-thick sources in general 
below. 

For the Milky Way (Sgr A*) we used the literature result 
from Baganoff et al. (2001) during quiescence. The data we 
use are displayed in Table 3 along with other galaxies with 
dynamically measured black holes without measurements of 
Lx, Lr, or either. A summary of the X-ray analysis may be 
gleaned from Figure 3, which shows a histogram of values 
of Eddington fractions /em = Lx /Lem for all objects that 
resulted in an X-ray measurement. The distribution shows 
that while most are accreting at a small fraction of Edding- 



ton, there are still a wide range of values encompassed in the 
sample. 

4. ANALYSIS 

4.1. Fitting Method 

For our measurement of the relation between Mbh, Lr and 
Lx, we considered the form 

l0gLfl, 38 =Rq + ^„logM BH ,8 + ^logL x ,40, (1) 

where we have normalized to Lr = 10 38 ergs^L^s, Mbh = 

10 8 M : Mbh,8; an d Lx = 10 40 Lx,40 m order to minimize in- 
tercept errors. To find the multi-parameter relation, we mini- 
mized the following statistic 

-2 _ ( R i - R Q- ~ ZpcXj) 2 

■ a" ■ + t z a" • + t z c ■ 

where R = logL^ 3g, p = logMBH.8, X = logLx,40, an d the 
sum is over each galaxy. The a terms are scatter terms that re- 
flect deviation from the plane due to intrinsic scatter and mea- 
surement errors. This statistic is the same statistic used by 
Merloni et al. (2003). We considered two cases. For the first, 
we assume that the intrinsic scatter is dominant and isotropic 
and thus use a total scatter projected in to the R direction: 
°o = a ri + %m a mi + ^°T.r To determine Co, we use a trial 
value of Co and increase the value until the reduced % 2 is 
unity after fitting with the new value. For the second, we use 
the measurement errors in Mbh an d Lx, assumed to be nor- 
mally distributed in logarithmic space, for o„, and a x respec- 
tively. The measurement errors in Lr are likely the smallest, 
and thus intrinsic scatter is likely to dominate. Here we as- 
sume a r = Oq. In this final case, our fit method is no longer 
symmetric, but it includes measurement errors and does not 
assume that the intrinsic scatter is isotropic. Both methods 
give nearly identical results, and we report only results from 
the latter method, which includes measurement errors. The er- 
rors on fit parameters come from the formal covariance matrix 
of the fit. 

4.2. Fundamental Plane Slopes 
Our best-fit relation for the fundamental plane is 

R = -0.34 ±0.24 
£ m = 0.78 ±0.27 

0.67 ±0.12. (3) 

The scatter we find in the Lr direction is Oq — 1 .00 dex, equiv- 
alent to 0.70 dex normal to the plane. These results are con- 
sistent with the findings of Merloni et al. (2003) and of Falcke 
et al. (2004). We plot several views of the fundamental plane 
in Figure 4 and the edge-on view in Figure 5. It is also inter- 
esting to note that for a fixed value of Mbh our relation finds 
Lr ex; L x 61 , consistent with the findings of Gallo et al. (2003). 

4.3. Mbh os the Dependent Variable 

We are using black hole masses that have been measured di- 
rectly. This approach allows us to use Lr and Lx as predictor 
variables for Mbh- We perform a multivariate linear regres- 
sion on Lr and Lxby assuming a form 

logM BH ,8 = W + c,-logL ff , 38 + c x \ogL X A0 (4) 
and minimizing 
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FIG. 1 . — Example Chandra spectra with best-fit models. The models have been folded through the instrument response. The horizontal error bars show the 
binning used for the fits. These four galaxies were chosen to show a variety of different models used to fit the data. All spectra included Galactic absorption and a 
power-law component. NGC 3031, NGC 4151, and NGC 4594 included intrinsic absorption; IC 4296 and NGC 4151 included an APEC model; and NGC 4151 
included a Gaussian component to model the Fe line. 




38 40 
\og(L x I erg s" 



46 



) [Literature] 



FIG. 2. — Comparison of results of X-ray analysis in this work to results 
from the literature. All values have been scaled to our adopted distances. 
Squares indicate straight-forward comparisons. Diamonds indicate that we 
have converted the literature result to an unabsorbed 2-10 keV luminosity 
using the published spectral fit and absorption. 
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FIG. 3. — Histogram of Eddington fractions defined as f^d = LxjL^&&. The 
contribution to the histogram from Seyfert galaxies is colored red, from other 
SMBH sources is colored blue, and from stellar-mass sources considered in 
section 5.3 is colored gray. The galaxy with the smallest /em is Sgr A*. A 
wide range of values are present in the sample even if most are found between 
/Edd — 10~ 9 and 10~ 6 . As expected, galaxies classified as Seyferts are, on 
average, emitting at a higher fraction of Eddington than other sources, and 
the stellar-mass sources are emitting at a higher fraction still. 



where o,„ is the measurement error in Mbh and Oo is an in- 
trinsic scatter term in the log (Mbh) direction. As before, the 
intrinsic scatter term is increased until the resulting best fit 



The Fundamental Plane with Dynamical Masses 



5 





FIG. 4. — Four views of the fundamental plane. Data are as described in Sections 2 and 3. Red points are galaxies classified as Seyferts. Blue points are 
LLAGNs and LINER galaxies. The varying views clearly show that as a whole the points lie on a plane in the dimensions shown. It is especially clear in the 
top-right panel that the LLAGN/LINER subsample appear to lie on a one-dimensional manifold. 



gives % 2 = 1 . We find a best-fit relation of 

fi = 0.19 ±0.19 
c r = 0.48 ±0.16 

c x = -0.24 ±0.15, (6) 

with an intrinsic scatter of Go = 0.77 dex in the mass direction. 
The intrinsic scatter is larger than other scaling relations (e.g., 
Oo = 0.44 ± 0.06 for the M-o relation and o = 0.38 ± 0.09 
for the M-L relation; Giiltekin et al. 2009b). We plot projec- 
tions of fit in the left panel of Figure 6. 

5. DISCUSSION 



5.1. Using A Black Hole 's Luminosity to Estimate Its Mass 

By using a sample of galaxies that have directly measured 
black hole masses, we are able to investigate the correlation 
between X-ray and radio luminosity and black hole mass. The 
measure of any correlation's worth as a predictor is the scatter, 
and we consider the scatter here. The scatter in the full relation 
is considerable (0.77 dex = 5.9), but it is only a factor of a cou- 
ple larger than other scaling relations used to estimate black 
hole mass. For example the M-a and M-L relations that relate 
Mbh and host galaxy velocity dispersion and bulge luminosity 
have intrinsic scatters of 0.44 dex = 2.75 and 0.38 dex = 2.4, 
respectively (Giiltekin et al. 2009b). 

It is worth noting that if we restrict the sample to just black 
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FIG. 5. — The fundamental plane relation. This figure shows the edge-on 
view of our best-fit relation: !;,„ = 0.78 and ^ v = 0.67. Error bars on the x-axis 
are calculated as aj = £%,oL , + 5v c 5 / ■ This view is primarily for comparison 
with Merloni et al. (2003) and with Falcke et al. (2004). Red circles are 
Seyferts. Blue circles are LINERs and unclassivied LLAGN. 

holes with mass M B h > 3 x 10 7 or M B h > 10 s M s , the intrin- 
sic scatter drops to Go = 0.45 or 0.41, respectively. There are 
several possible interpretations for the decreased scatter when 
restricting the sample by mass. One possibility is that the re- 
quirement of detection in both radio and X-rays translates to 
a requirement of high Eddington fraction for low-mass black 
holes at a fixed distance. The mean values of /gdd for the 
whole sample, for the sample with Mbh > 3 x 1O 7 M , and for 
the sample with Mbh > 10 8 M ... are approximately 6 x 10~ 5 , 
6 x 10~ 6 , and 3 x 10~ 6 , respectively. It is possible that when 
sources accrete at a higher rate, the fundamental plane relation 
may no longer apply. 

Another possible explanation for the smaller scatter in the 
high-mass sample is that the low-scatter trend is real, and that 
the scatter estimated from the entire sample is skewed by a 
few data points. The most obvious outliers from the left panel 
of Figure 6 are Circinus and NGC 1068. If these two are elim- 
inated, the scatter becomes Go = 0.50 dex. The derived intrin- 
sic luminosities of these sources may be difficult to determine 
because of obscuration. In these sources we have a poor view 
of the central engine and are seeing reflected, rather than direct 
X-ray emission (Matt et al. 1996; Antonucci & Miller 1985). 
If the intrinsic X-ray luminosity of these sources is higher, 
then they would lie closer to the best-fit plane than they do 
now. 

AGN classification for each galaxy of the sample is listed 
in Table 3. The distinction between Seyferts and LINERs is 
judged from the line ratios with the usual diagnostic and divi- 
sion set so that Seyferts have [OIII]A.5007/Hp > 3.0 (Veilleux 
& Osterbrock 1987) as a measurement of the level of nuclear 
ionization, though there is no obvious transition between the 
two classes (Ho et al. 2003). The physical difference between 
LINERs and Seyferts may be that the LINERs lack a "big 
blue bump" and produce a larger partially ionized zone. The 
transition in spectral energy distribution from a Seyfert to a 
LINER may happen at low /gdd (Ho 2008). The distinction 
between Seyfert types is determined by the ratio of broad-line 
and narrow-line emission. LLAGN are defined by having an 
Ha luminosity smaller than L(Ha) < 10 40 ergs -1 (Ho et al. 
1997). The difference between Seyfert types is understood to 
be due to differing viewing angles with respect to an obscuring 
dusty torus that surrounds the broad line region (with type 1 
unobscured and type 2 completely obscured). For a review of 



the observational differences among the different classes and 
the current physical explanations for the differences see the 
review by Ho (2008). 

We may give special consideration to all non-Seyfert AGNs 
in our sample. Since all Seyferts in our sample are at least 
partially obscured, obscuration is one potential issue that is 
addressed. Obscuration will naturally lead to an underesti- 
mate in X-ray luminosity. We minimize this by fitting for 
the absorption across the 0.5-10 keV band. Since the softer 
photons are more readily absorbed, the shape of the spectrum 
gives an indication of the level of absorption. We also use the 
hard X-ray flux, which is least affected by absorption, for our 
X-ray luminosity. Nevertheless, the most heavily obscured 
sources may be intrinsically brighter than our fits indicate. 
We attempt to isolate this issue below by removing Compton- 
thick sources. In addition to obscuration, as mentioned above, 
Seyferts also accrete at higher fractions of Eddington and may 
accrete in a mode different from LINERs. In addition, since 
Seyferts are thought to be dominated by thermal output, their 
radio luminosities may be poor probes of the power in out- 
flows and thus not belong on the relation considered here. 
Thus, there is a physical motivation to separate them from the 
rest of the sample. 

When we only use the 8 LINER and unclassified LLAGN 
sources, our fit becomes 

/j = 0.70±0.11 
c r = 0.48 ±0.14 

c x = -0.07 ±0.12, (7) 

with a scatter of Go = 0.25, substantially smaller than other 
intrinsic scatter measurements found for this relation and ac- 
tually smaller than the scatter in the M-G and M-L relations. 
Kording et al. (2006) similarly found a substantially reduced 
scatter in fundamental plane fits to a sample of only stellar- 
mass black holes, Sgr A*, and LLAGNs. The fit we find is 
significantly different from the other fits, notably that it is con- 
sistent with no dependence on X-ray luminosity (c x = 0). This 
is at odds with the findings of Ho (2002), who found no de- 
pendence of black hole mass on radio luminosity. The data 
do appear to lie on a one-dimensional manifold in the three- 
dimensional space considered, but with only 8 data points, the 
data set is substantially smaller than that of Ho (2002), who 
also used direct, primary mass measurements in addition to di- 
rect, secondary mass measurements (i.e., reverberation map- 
ping). 

If obscuration, rather than accretion rate or accretion 
mode, is the underlying reason for the smaller scatter in the 
LINER/LLAGN sample, then we should see similar results 
when omitting sources that are Compton thick («# >Gj — 
1 .5 x 10 24 cm -2 ). Compton-thick sources will be heavily ob- 
scured and the intrinsic luminosities may be much higher than 
the observed flux would imply (Levenson et al. 2002, 2006). If 
we conservatively omit the sources from Table 2 intrinsic ab- 
sorption larger than 10 24 cm -2 (NGC 3031, NGC 4374, and 
NGC 6251) as well as the sources determined to be Compton 
thick from Fe Ka modeling (Circinus and NGC 1068; Leven- 
son et al. 2002, 2006), we obtain 

/jo = 0.40 ±0.16 
c r = 0.46 ±0.13 

c r = -0.14±0.12, (8) 

with a scatter of Go = 0.53. This result is consistent with the 
Seyfertless sample at about the la level, though with a larger 
scatter. 




FIG. 6. — Best fit linear regression of Mbh on Lr and L\ for (left) all galaxies and for (right) LLAGN and LINER galaxies only. The relation on the right is 
considerably tighter but may be affected by the small number of sources. Red circles are Seyferts. Blue circles are LINERs and unclassivied LLAGN. 



5.2. SgrA* 

Sgr A*, the central black hole in the Galaxy, is a unique 
source in many ways. Its extremely low accretion rate 
(Lx/L^dd ~ 4 x 10 -12 ) is two orders of magnitude below the 
next lowest in our sample. An analog to Sgr A* could not be 
observed outside of the local group. 

When using only the two nearby super-massive black holes 
with extremely well-determined mass and distance (Sgr A* 
and NGC 4258) and the X-ray binary in which the correla- 
tion extends over several orders of magnitude (GX 339—4), 
the best-fit fundamental plane relation changes so that Sgr A* 
is under-luminous in X-rays during quiescence by at least 2 
orders of magnitude (Markoff 2005). Such a break from the 
correlation (also seen in some X-ray binaries as they rise out 
of quiescence Coriat et al. 2009) may indicate that, during qui- 
escence at least, Sgr A* is accreting in a different mode than 
the correlation sources. If such an extremely low accretion 
rate is in a different category from the rest of the objects, then 
it makes sense to exclude Sgr A* from the sample, in which 
case our best fit becomes 

R = -0.46 ±0.30 
$ m = 0.84 ±0.29 

E, x = 0.80 ±0.23, (9) 

with an intrinsic scatter of Oo = 1.02, which is not a signifi- 
cantly different fit. 

5.3. Stellar-mass Sources 

Our initial sample includes only the supermassive black 
holes in galactic centers. There are, however, several Galactic 
stellar-mass black holes with dynamically measured masses. 
If accretion onto black holes is driven by the same physi- 
cal processes at all mass scales, then the stellar-mass sources 
should obey the same relation, which is what Merloni et al. 
(2003) and Falcke et al. (2004) found. So while our focus has 
been on super-massive black holes, we may revisit our calcu- 
lations with the sample of stellar-mass black holes given in 
Table 1. This sample was selected from stellar-mass black 
holes with dynamically determined masses with simultane- 
ous X-ray and radio data. In addition to the sources listed, 
there were two stellar mass black holes that had adequate data 
(4U 1543-475 and GRO J1655-40) but whose jets may not 
be in a steady state and thus skewing the relation. 

The stellar-mass systems, with the possible exception of 
GRS 1915, are in the low/hard state, which is characterized 
by a hard X-ray photon index (1.4<T<2.1),a small ratio of 



unabsorbed disk flux to total unabsorbed flux (/ < 0.2; Remil- 
lard & McClintock 2006) and is usually seen at low Eddington 
rates. This state is also typically associated with a steady ra- 
dio jet whereas jets in the high/soft state are quenched (Fender 
2001). By requiring radio emission, we essentially require 
a low/hard state. If such a state can be extended to SMBH 
sources, it would naturally compare with the similarly low 
Eddington rates in LLAGNs in which jet emission is more 
prominent compared to Seyferts. The mapping of X-ray bi- 
nary states to accreting SMBHs is complicated by the fact that 
no comparable transitions are seen in SMBHs. 

These three accreting black holes have masses measured 
from period measurements of the donor star's orbit. The mass 
of the donor star is estimated based on spectral type, and the 
inclination of the orbit for systems such as these is generally 
derived from modeling the star's change in flux, assumed to 
be from the change in viewing angle of a tear-drop-shaped 
object (ellipsoidal modulation). For two of the three stellar- 
mass sources we are using, however, the inclination is con- 
strained by other means. For GRS 1915+105 the inclination is 
constrained from the apparent superluminal motion of ejected 
jet material that is assumed to be perpendicular to the orbital 
plane based on the lack of observed precession (Mirabel & 
Rodriguez 1994; Greiner et al. 2001b). For Cyg X-l, the in- 
clination has been estimated in several ways, including UV 
line modeling and X-ray polarization (Ninkov et al. 1987, and 
references therein). 

The luminosity data from each source is simultaneous, 
which is important for these highly variable sources. For two 
of the sources, we use two sets of simultaneous observations. 
Using more than one observation of a particular source in the 
fit over-weights that source and will skew the fit if it is atypi- 
cal. Under the assumption that each source belongs in the fit 
in all of the epochs used, however, they provide valuable extra 
information of possible accretion states in the same relation. 

The results of our fundamental plane fits become: 

R = -0.33 ±0.21 
$ m = 0.82 ±0.08 

^= 0.62 ±0.10, (10) 

with an intrinsic scatter of Oo = 0.88. The uncertainties in 
slopes have decreased because of the increased range in the 
values present, especially for t, m . It is interesting to note that 
while the best-fit parameters do not significantly change from 
our fits to central black holes, the intrinsic scatter does. This 
decrease can be attributed to the fact that these sources lie 
closer to the plane. It is also worth noting that the fits do not 
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TABLE 1 
Stellar-Mass Black Hole Data 



Name 


D 


log(M BH ) 


log (I*) 


log(Ljr) 


Refs. 


GRS 1915+105 


11 


1.15±0.13 


30.64 a 


38.06±0.06 


1,2,3,3 




11 


1.15±0.13 


30.90 a 


38.69±0.06 


1,2,3,3 


V404 Cyg 


3 


1.08±0.07 


28.30 


33.07 ±0.24 


4,4,5,6 


Cygnus X-l 


2.5 


1.00±0.24 b 


29.91 c 


36.71±0.18 d 


7,8,9 




2.5 


1.00±0.24 b 


29.84 c 


36.77±0.18 d 


7,8,9 



of scatter in the relation may be an indication that they are 
diverging away from the fundamental plane relation. To bet- 
ter understand the differences between Seyfert galaxies and 
the other sources, a future theoretical work will consider just 
these types of sources, including physical modeling of the data 
sets presented here. 

6. CONCLUSIONS 

In this paper we analyze the relationship among X-ray lumi- 
References. — (1) Fender et al. 1999; (2) Greiner et al. 2001a; (3) Muno et absity, radio luminosity, and the mass of a black hole. Distinct 

2001; (4) Shahbaz et al. 1994; (5) Gallo et al. 2005a; (6) Bradley et al. 2007; (g^ m p rev ious Studies of this relationship, we use only black 
Bregman et al. 1973; (8) Herrero et al. 1995; (9) Stirling et al. 2001. . . . . i_ j • .. j n r 

b 6 hole masses that have been dynamically measured. Because of 

NOT t E ; ~~ of U f TT s black hole data , used in Sf l cti °. n 53 - Distances are given ^ re i at ively small distances to the objects in this sample, we 

units ot kpc. Black hole masses are in solar units. Radio and X-ray luminosities are . . . , , . . . „ a J ,. , r 

' -l ah i . , , , . avoid potential biases arising from flux limited samples. Us- 

ui units ot erg s . All values are scaled to the distances given. I he sources were .in r r , , 

low/hard state for the epochs listed with the exception of GRS 1915, which may be iAiag the most recent compilation of black hole masses, we ana- 
plateau state (Muno et al. 2001). The numbers in the reference column give the numliyzed archival Chandra data to get nuclear X-ray luminosities 

of the original reference for the distance, mass, radio luminosity, and X-ray luminosijyj jjjg J7 — 2—10 keV band We combined this with V = 5 GHz 
respectively. X-ray luminosities have been converted to the E = 2-10 keV band. ^ luminosities found in the litera ture and fit a relation of 
a Interpolated from v = 2.25 GHz and v = 8.3 GHz data. . form L R 38 = R + £ m logM BH 8 + UogZ* 40 to find 

Mass uncertainty was estimated from the range ot values found in the literature (Mc- ' * ' 

Clintock & Remillard 2006). n n _ . , ~ „ ■ 

Kq = — U. 34 ± U.z4 



c Extrapolated from V = 8.4 GHz assuming constant vF v . 
d Data come from Rossi X-ray Timing Explorer (RXTE) All-Sky Monitor (ASM) as 
suming a standard spectral form. 



change even though two of the stellar-mass sources are accret- 
ing at a much higher fraction of Eddington than the supermas- 
sive sources. GRS 1915+105 is accreting at/gdd ~ 0.06 to 0.3, 
and Cygnus X-l at / Edd w 0.004 to 0.005, whereas all of the 
supermassive sources are accreting at /edd < 0.001 (Fig. 3). 

It should be noted that there are different systematic er- 
rors in the stellar-mass and central black holes. The mass 
measurements are from completely different methods. The 
X-ray extragalactic sources may be contaminated from point 
sources and may be more heavily obscured than the stellar- 
mass sources. The extragalactic sources may also be contam- 
inated by supernova remnants along the line of sight, though 
this can be mitigated by going to higher frequencies. Stel- 
lar mass uncertainties are dominated by uncertainties in dis- 
tance, inclination, and light from accretion (see Reynolds et al. 
2008). 

5.4. Future Work 

In this paper, we have only included the 18 black holes with 
measured masses, radio fluxes, and X-ray fluxes. This sam- 
ple makes up slightly more than one third of the entire sample 
of black holes with measured masses. There are 1 1 without 
nuclear radio data or only with upper limits on one or more 
of these quantities. There are a further 16 sources with no 
Chandra X-ray fluxes measured because either there are no 
Chandra data or merely insufficient data. Many of the sources 
have masses M < 1O 8 M . By completing the sample of M-a 
black holes with further X-ray and radio observations, the in- 
creased number of data points should be especially helpful in 
determining whether the large scatter at the low-mass end and 
the small scatter at the high-mass end are actual differences or 
just artifacts of a few outliers. 

Another place for future work is in understanding the ap- 
parent special place that Seyfert galaxies occupy in the fun- 
damental plane. If one were to naively assign accretion 
states used for stellar-mass black holes to Seyfert galaxies, 
they would be considered in the thermally dominant/high-soft 
state. For stellar-mass black holes in this state, jets are not 
measured. That the Seyfert galaxies are an apparent source 



0.78 ±0.27 
0.67±0.12 



(11) 



with a scatter of a = 1.00 in the logL^ direction, consistent 
with previous work. We also fit a relation to be used as an 
estimation for black hole mass based on observations of Lx 
and Lftof the form: 



logM BH ,8 = A<0 + C,-l0gLfl.3g + C A -logLx,40, 



(12) 



finding 



0.19±0.19 
0.48 ±0.16 
-0.24±0.15, 



(13) 



with an intrinsic scatter of Oo = 0.77 in the logMeH direc- 
tion. This intrinsic scatter is larger than other scaling relations 
involving Mbh, but decreases considerably when only using 
the most massive black holes or when eliminating obscured 
central engines from the sample. Both of these issues require 
further investigation and could be answered by completing the 
sample with more Chandra observations. 
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TABLE 2 
Summary of Chandra Spectral Fits 
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Note. — 


Results from X-ray spectral analysis 


. First column gives galaxy name. The second column 


gives Chandra observation identification number. 


The third column lists exposure time in units of ks. Fourth column lists % jv 



where V is the number of degrees of freedom. If the fit used C-stat statistics instead of % 2 statistics, then the third column is left blank. Best-fit parameters with la errors for each. A blank entry in a given column indicates that the given 
component was not part of the spectral model used. Galaxies with superscript "a" were only able to constrain an upper limit to the flux. The model for Circinus also included a pileup model. 



Giiltekin et al. 



TABLE 3 
Black Hole Data 



fulfil n y v 


AGN Class 


n/Mnr 


1 OCT I Mr,,, 1 


\CIP (J v\ 




Ref. 


Circinus * 


S2 a 


4.0 


6.23±0.088 


37.73 


41.48±0.034 


1,2 


IC 1459 * 


S3 


30.9 


9.44±0.196 


39.76 


40.86±0.014 


3.4 


IC 4296 * 




54.4 


9.13±0.065 


38.59 


41.31±0.044 


5.6 


Sgr A* * 




0.008 


6.61±0.064 


32.48 


33.33±0.068 


7.S 


NGC 0221 




0.9 


6.49±0.088 




36.17±0.059 


9 


NGC 0224 


S3" 


0.8 


8.17±0.161 


32.14 


< 36.00 


10,11 


NGC 0821 




25.5 


7.63±0.157 




38.44±0.640 


12 


NGC 1023 




12.1 


7.66±0.044 




38.80±0.066 


13 


NGC 1068 * 


S2 J 


15.4 


6.93±0.016 


39.18 


39.54±0.024 


14,15 


NGC 1300 




20.1 


7.85±0.289 






16 


NGC 1399 




21.1 


8.71 ±0.060 




< 38.64 


17 


NGC 2748 




24.9 


7.67±0.497 






16 


NGC 2778 




24.2 


7.21±0.320 






12 


NGC 2787 * 


S3b 


7.9 


7.64±0.05() 


36.52 


38.7O±0.O59 


18,19 


NGC 3031 * 


S1.8 a 


4.1 


7.90±0.087 


36.97 


4().84±0.097 


20,21 


NGC 3115 




10.2 


8.98±0.182 




38.04±0.312 


22 


NGC 3227 * 


S1.5 


17.0 


7.18±0.228 


37.72 


41.55±0.046 


23,21 


NGC 3245 * 


S3 11 


22.1 


8.35±0.106 


36.98 


39.28±0.420 


23,24 


NGC 3377 




11.7 


8.06±0.163 




38.00±0.322 


12 


NGC 3379 * 


S3" 


11.7 


8.09±0.250 


35.81 


38.17±0.205 


25,26 


NGC 3384 




11.7 


7.25±0.042 




< 38.55 


12 


NGC 3585 




21.2 


8.53±0.122 




38.98±0.161 


27 


NGC 3607 


S2 


19.9 


8.08±0.153 




< 38.60 


27 


NGC 3608 


S3' 1 


23.0 


8.32±0.173 




< 38.79 


12 


NGC 3998 * 


S3b 


14.9 


8.37±0.431 


38.03 


41.44±0.007 


28,29 


NGC 4026 




15.6 


8.33±0.109 




< 38.53 


27 


NGC 4258 * 


S2 


7.2 


7.58±0.001 


36.03 


40.83±0.096 


30,21 


NGC 4261 


S3h 


33.4 


8.74±0.090 


39.32 


< 40.92 


31,29 


NGC 4291 




25.0 


8.51 ±0.344 






12 


NGC 4342 




18.0 


8.56±0.185 




39.13±0.151 


32 


NGC 4374 * 


S2 


17.0 


9.18±0.231 


38.77 


39.42± 1.503 


33,34 


NGC 4459 


S3' 1 


17.0 


7.87±0.084 


36.13 


< 38.97 


18,24 


NGC 4473 




17.0 


8.11±0.348 




< 38.50 


12 


NGC 4486 * 


S3 


17.0 


9.56±0.126 


39.83 


40.46±0.015 


35,36 


NGC 4486A 




17.0 


7.13±0.146 




< 38.96 


37 


NGC 4564 




17.0 


7.84±0.045 




< 38.79 


12 


NGC 4594 * 


S1.9 


10.3 


8.76±0.413 


37.89 


40.19±0.307 


38,39 


NGC 4596 


S3 a 


18.0 


7.92±0.162 




< 38.72 


18 


NGC 4649 




16.5 


9.33±0.117 


37.45 


< 38.95 


12,40 


NGC 4697 




12.4 


8.29±0.038 




38.25±0.745 


12 


NGC 5077 


S3b 


44.9 


8.90±0.221 






12 


NGC 5128 * 


S2? 


4.4 


8.48±0.044 


39.85 


40.22±0.085 


41,42 


NGC 5576 




27.1 


8.26±0.088 






27 


NGC 5845 




28.7 


8.46±0.223 




39.07±0.722 


12 


NGC 6251 * 


S2 


106.0 


8.78±0.151 


41.01 


42.50±0.207 


43,44 


NGC 7052 




70.9 


8.60±0.223 


39.43 


< 40.69 


45,46 


NGC 7457 




14.0 


6.61±0.170 




< 38.28 


12 


NGC 7582 * 


S2 a 


22.3 


7.74±0.104 


38.55 


41.69±0.208 


47,48 


PGC 49940 




157.5 


9.59±0.056 






5 


Cygnus A 


S1.9 


257.1 


6.43±0.126 


41.54 


44.23±0.088 


49,6 


NGC 4151 


S1.5 


13.9 


7.65±0.048 


38.20 


41.69±0.074 


50,21 


NGC 4303 


S2 


17.9 


6.65±0.349 


38.46 


38.89±0.124 


51,52 


NGC 4742 




16.4 


7.18±0.151 






53 


NGC 4945 


S 


3.7 


6.15±0.184 


38.17 


37.80±0.921 


54,55 


NGC 5252 


S2 


103.7 


9.00±0.341 


39.05 


43.20±0.017 


56,57 
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TABLE 3 
Black Hole Data 



References. — (1) Greenhill et al. 2003, (2) Turner & Ho 1983, (3) Cappellari et al. 2002, (4) Sadler et al. 1989, (5) Dalla Bonta et al. 2008, (6) Sambruna et al. 1999, (7) Ghez 
et al. 2008 and Gillessen et al. 2008, (8) Ekers et al. 1983, (9) Verolme et al. 2002, (10) Bender et al. 2005, (11) Crane et al. 1992, (12) Gebhardt et al. 2003, (13) Bower et al. 2001, 
(14) Lodato & Berlin 2003, (15) Ulvestad & Wilson 1984, (16) Atkinson et al. 2005, (17) Gebhardt et al. 2007, (18) Sarzi et al. 2001, (19) Heckman et al. 1980, (20) Devereux et al. 
2003, (21) Ho & Ulvestad 2001, (22) Emsellem et al. 1999, (23) Barth et al. 2001, (24) Wrobel & Heeschen 1991, (25) Gebhardt et al. 2000, (26) Fabbiano et al. 1989, (27) Giiltekin 
et al. 2009a, (28) de Francesco et al. 2006, (29) Wrobel & Heeschen 1984, (30) Herrnstein et al. 2005, (31) Ferrarese et al. 1996, (32) Cretton & van den Bosch 1999, (33) Bower et al. 
1998, (34) Jenkins et al. 1977, (35) Macchetto et al. 1997, (36) Biretta et al. 1991, (37) Nowak et al. 2007, (38) Kormendy 1988, (39) Hummel et al. 1984, (40) Spencer & Junor 1986, 
(41) Silge et al. 2005, (42) Wright et al. 1994, (43) Ferrarese & Ford 1999, (44) Jones et al. 1986, (45) van der Marel & van den Bosch 1998, (46) Morganti et al. 1987, (47) Wold et al. 
2006, (48) Gregory et al. 1994, (49) Tadhunter et al. 2003, (50) Onken et al. 2007, (51) Pastorini et al. 2007, (52) Gregory & Condon 1991, (53) listed as in preparation in Tremaine 
et al. 2002 but never published, (54) Greenhill et al. 1997, (55) Elmouttie et al. 1997, (56) Capetti et al. 2005, (57) Polletta et al. 1996. 

NOTE. — This table lists all galaxies with dynamically measured black hole masses. Sources with an asterisk after their name are those used in this paper. The second column gives 
AGN classification from Veron-Cetty & Veron (2006) unless it has a superscript "a," in which case it comes from NED. "SI" indicates type 1 (unobscured) Seyfert, "S2" indicates 
type 2 (obscured) Seyfert, "Sl.X" indicates transitional or intermediate Seyfert, "S3" indicates type 3 Seyfert or LINER galaxy, and "?" indicates that NGC 5128 is a questionable 
BL Lac object. Note that NGC 3227 and NGC 4151 are classified as type 1.5 but both have reverberation mapping masses (Onken et al. 2003; Bentz et al. 2006) and thus have visible 
broad-line regions. Beyond these two galaxies, none of the galaxies is obviously a Seyfert 1, though NGC 1068 and NGC 7582 are classified as such by Veron-Cetty & Veron (2006). 
We classify them according to their NED classifications as Seyfert 2 galaxies. NGC 1068 is a Seyfert 2 galaxy and only shows broad Balmer lines in polarized light, indicating that the 
light has been scattered and thus coming from behind an obscured source (Antonucci & Miller 1985). NGC 7582 is a classical Seyfert 2 galaxy that developed broad emission lines 
for a short period of time in 1998 July (Aretxaga et al. 1999). The change to a Seyfert 1 spectrum may be explained by a stellar disruption event, a change in the obscuring medium, or 
a type Iln supernova (Aretxaga et al. 1999, see also Veron-Cetty & Veron 2000). It has also been suggested, based on its X-ray spectrum, that NGC 7582 is an obscured narrow-line 
Seyfert 1 (Dewangan & Griffiths 2005). For our purposes, we classify this source as a Seyfert 2. The paucity of true Seyfert 1 galaxies in our sample is not surprising since such bright 
central engines would compromise dynamical black hole mass measurements. The third column gives distance to the galaxy in units of Mpc, which is used to scale all data. The 
fourth column lists logarithmic black hole mass per unit solar mass as compiled by Giiltekin et al. (2009b). The fifth column gives logarithmic radio luminosity in units of ergs -1 . The 
radio data come from the compilation of Ho (2002) with the following exceptions: IC 4296, NGC 5128, NGC 7582, NGC 4303, and NGC 5252. The sixth column gives logarithmic 
X-ray luminosity in units of ergs -1 , which come from this work except for Sgr A* (Baganoff et al. 2001) and the upper limit on NGC 0224 (Garcia et al. 2005). We leave the column 
blank if there are no archival data available. The final column lists original references for the Mbh measurement and radio luminosity, if present. The bottom portion of the table gives 
data for when the black hole mass may be wrong. For this paper we only use data from galaxies that are in the top portion and that have both radio and X-ray detections. Blank entries 
indicate that there were no archival data available and may be followed up with more observations. 



